Pairing Fluctuations and Anomalous Transport Above the BCS-BEC Crossover 
in the Two Dimensional Attractive Hubbard Model 



O 

(N 



in 



I 



I 

O 

o 



> 
in 



in 
o 



X 



Sabyasachi Tarat and Pinaki Majumdar 

Harish- Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211019, India 

(Dated: 5 May 2011) 

A Fermi liquid with weak attractive interaction undergoes a BCS transition to a superconductor 
with reducing temperature. With increasing interaction strength, the thermal transition is pro- 
gressively modified as the high temperature 'metallic' phase develops a pseudogap due to pairing 
fluctuations and the resistivity above Tc shows insulating behaviour. The crossover to insulating 
character occurs much before the system can be considered to be in the BEC regime of preformed 
fermion pairs. We use a new Monte Carlo tool to map out the BCS-BEC crossover in the attractive 
Hubbard model on large two dimensional lattices and explicitly compute the resistivity to demon- 
strate how the metal to superconductor (MS) thermal transition at weak coupling crosses over to 
an insulator to superconductor (IS) transition at intermediate coupling. Our high resolution access 
to the single particle and optical spectrum at finite temperature allows us to completely describe 
the transport crossover in this longstanding problem. 



The BCS-BEC crossover in attractive fermion systems 
has been a topic of interest [ij for several decades. With 
increasing interaction, the ground state of a weak cou- 
pling 'BCS superconductor', with pair size ^ much larger 
than the interparticle separation ~ fc^^ (where fc^? is the 
Fermi wavevector) evolves smoothly [^tSI into a 'Bose- 
Einstein condensate' (BEC) of preformed fermion pairs 
with ^ < fc^^. The 'high temperature' normal state 
changes from a conventional Fermi liquid at weak cou- 
pling to a gapped phase at strong coupling. While the 
pairing gap increases with coupling strength, the super- 
conducting Tc in lattice models reaches a maximum at in- 
termediate coupling and falls thereafter. A striking con- 
sequence of the separation of pairing and superconduct- 
ing scales is the emergence of a (pseudo)gapped normal 
phase, with preformed fermion pairs but no superconduc- 
tivity due to strong phase fluctuations. 

The early work of Leggett Q and Nozieres and 
Schmitt-Rink provided the intuitive basis for under- 
standing this problem. It has since been followed up 
by extensive quantum Monte Carlo (QMC) work [71- 
powerful semi-analytic schemes |12h14|- and most 
recently dynamical mean field theory (DMFT) [11, [3. 
The efforts have established the presence of a pseudogap 
[^ in the single particle spectrum beyond moderate cou- 
pling and temperature T > T^, and also a gap in the spin 
excitation spectrum [T^j- While these indicate a break- 
down of the Fermi liquid picture, the crucial transport 
properties in the normal state remain obscure. 

For example, given the 'gapped' normal state at strong 
coupling, is it insulating? If so, down to what coupling 
does this extend? There seem to be various options: 
(i) the insulating state arises only when the single parti- 
cle density of states (DOS) at T > Tc has a 'hard gap', 
or (ii) metallic conduction survives even in the hard gap 
state due to transport by 'bosonic' carriers, or (iii) the 
normal state becomes insulating, due to strong 'pairing 
disorder', even without a hard gap in the spectrum. QMC 



calculations, which set the benchmark in the field, unfor- 
tunately do not have access to real frequency information 
or the system size that can resolve this issue. 

We use a new Monte Carlo method, involving a 
Hubbard-Stratonovich (HS) decomposition [17,] of the at- 
tractive Hubbard model in terms of pairing fields, in two 
dimensions (2D). We ignore the time fluctuations of the 
(bosonic) pairing field treating it as classical [l^, but 
fully retain the spatial fluctuations of its amplitude and 
phase. Our principal results are the following: (1). We 
benchmark the Tc obtained by our method with the most 
recent QMC results, demonstrating the accuracy of our 
method at temperatures of interest. (2). We obtain the 
resistivity p(T) over the entire coupling range, and ob- 
serve that the high temperature phase goes insulating at 
U/t > 3, in the nominally 'BCS regime' and much before 
the occurence of preformed pairs. (3). The single par- 
ticle DOS reproduces features previously inferred from 
QMC, but the optical spectrum behaves very differently 
from the single particle DOS, 'filling up' at a much lower 
temperature at strong coupling. 

We study the attractive Hubbard model in 2D. 
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t denotes the nearest neighbour tunneling amplitude, 
[/ > is the strength of onsite attraction, and ji the 
chemical potential. We will study the range 1 < U /t < 
10, going across the BCS-BEC crossover, and set ^ so 
that the particle density remains at n w 0.9. 

The model is known to have a superconducting (S) 
ground state for all n ^ 1, while at n = 1 there is the 
coexistence of superconducting and density wave (DW) 
correlations in the ground state. For n ^ 1 the ground 
state evolves from a BCS state at J7/t ^ 1 to a BEC of 
'molecular pairs' at {//t ^ 1. The pairing amplitude and 
gap at T = can be reasonably accessed within mean 
field theory or simple variational wavefunctions. 
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Mean field theory, however, assumes that the elec- 
trons are subject to a spatially uniform self-consistent 
pairing amplitude ((c|t4i))- small U/t this vanishes 
when ksT ^ ie^*/*^, but at large U/t it vanishes only 
when ksT ^ U. The actual Tc at large U is con- 
trolled by phase correlation of the local order parame- 
ter, rather than finite pairing amplitude, and occurs at 
ksTc ~ f{n)t'^/U, where f{n) is a function of the den- 
sity. The wide temperature window, between the 'pair 
formation' scale fc^T/ ^ U and fc^Tc corresponds to equi- 
librium between unpaired fermions and hardcore bosons 
(paired fermions). QMC calculations have suggested a 
single particle pseudogap and a gap in the spin (NMR) 
excitation spectrum for U/t >; 4 and T > Tc. None of 
the calculations, however, seem to have addressed the 
simplest measurable property, i.e, normal state charge 
transport. This is probably due to 'analytic continuation' 
problems in QMC data or the severe finite size effects in 
exact diagonalisation based schemes. 

We use a strategy used earlier to access the supercon- 
ductor to insulator transition in the disordered attractive 
Hubbard model , and models of d — wave pairing [3| , 
augmented now by a Monte Carlo technique that allows 
access to system size upto ~ 40 x 40, much larger than 
the coherence length for the chosen U >2. We decouple 
the Hubbard term in the pairing channel by using the 
HS transformation and treat the HS field in the static 
approximation. This gives qualitatively correct answers 
at all n 7^ 1, and surprisingly accurate Tc values when 
compared to full QMC 20]. 

The static HS approach leads to the effective model: 
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where Hq = -tY.{^j)a ~ ^^■Y.^a^^'y and A, ^ 

Ai|e*^' is a complex scalar classical field. This model 
allows fiuctuations in both the amplitude and phase of 
the HS variable, and the fermions propagate typically in 
an inhomogeneous background defined by A^. 

To obtain the ground state, and in general configura- 
tions {|Aj|,0i} that follow the distribution P{|Aj|,0i} cx 
Q-l3H^ff^ we use the Metropolis algorithm to update 
the lAI and 9 variables. This involves solution of the 



Bogoliubov-de Gennes (BdG) equation [21| for each at- 
tempted update. For equilibriation we use a 'traveling 
cluster' algorithm [22j, diagonalising the BdG equation 
on a 8 X 8 cluster around the update site. Global proper- 
ties like pairing field correlation, DOS, etc, are computed 
via solution of the BdG equation on the full system. All 
results in this paper are for system size 32 x 32. 

If are the BdG eigenvalues in some equilibrium 
configuration a, the quasiparticle DOS is computed as 
N{uj) = S{uj — e")), where the angular brackets in- 
dicate averaging over a. The quasiparticle (QP) gap is 
the minimum of over all a, n at a given T. Similarly, 



the optical conductivity in an equilibrium background 
is formally cr(uj) = —uj~^Im{hxx(<l = 0,w)) where the 
current-current correlation function is defined by 
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The |n), \m) are multiparticle states of the system. We 
have suppressed the a labels above. We will discuss the 
simplification of this expression elsewhere. For the con- 
tinuous part of cr(w), i.e, excluding the superfiuid re- 
sponse, it leads to: 



ct{lo) = ^i^i(a,6) 
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where, now, the £0,6/3 > 0, etc, are single particle eigen- 
values of the BdG equations, the n(ea), etc., are Fermi 
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FIG. 1: Colour online: (a). The transition temperature T'c((7). 
(b). Superconductor- metal-insulator (S-M-I) character. The 
superconductor has pdc = 0, the M and I have pdc finite, with 
dpdc/dT > for M and dpdc/dT < for /. (c). 'Phases' 
according to single particle DOS near Fermi level: UG is 
ungapped (Fermi liquid), PG is pseudo-gapped, and G is 
gapped. By this criteria the superconductor and T > Tc 
insulator are similar, (d). Classification according to low fre- 
quency optical spectrum, i. e, the presence of the superfiuid 5 
function, and a continuous weight extending to cj = 0. Both 
PI and P2 are in the superfiuid region but also have low 
frequency optical spectral weight (except at T = 0). The 
distinction between PI and P2 is quantitative, P2 has much 
larger low u weight. P3 and P4 are both 'normal' i.e, with- 
out a superfiuid peak, and also have low uj weight. Again, 
the PS, P4 distinction is quantitative, the low uj weight is 
exponentially small in P4. 
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FIG. 2: Colour online: The resistivity p{T) for various inter- 
action strengths. There is a region with dp/dT < above Tc 
at all U, including the weak U systems. Note that the y axis 
spans 12 orders of magnitude. 



functions, and the F's are current matrix elements com- 
puted from the BdG eigenfunctions. We will call the two 
contributions above, aiicj) and (T2(w). The averaging of 
cr{uj) over equilibrium configurations of at a given T is 
implied. The dc resistivity p = ijJo~^ a{uj)du}, where 
Wo ^ O.lt. 

Phase diagram: Fig.l.(a) shows our result on Tc(U), 
at n — 0.9. We compute the thermally averaged pair- 
ing field correlation S'(q) -^{J^ij \^i\\^j\cos{9i - 

6'j)e^i('''~""^)) at q = {0,0}. This is like the 'ferro- 
magnetic' correlation between the A^, treating them as 
two dimensional moments. If the q — {0, 0} component, 
S{0), is 0(1) it implies that the pairing field has a non- 
zero spatial average and would in turn induce long range 
order in the thermal and quantum averaged correlation 
Mij = {{cl^cl^CjiCj^)). We locate the superconducting 
transition from the rise in S{0, T) as the system is cooled. 
The results are not reliable below U/t < 1, since ^ be- 
comes comparable to our system size, but compare very 
well with available QMC data for U/t ^ 2. 

Fig.l.(b) shows the basic classification of the U — T 
space in terms of metal (M), insulator (I) and super- 
conductor (S), the M and I regions determined from the 
slope dp/dT (see Fig. 2). The remarkable feature is that 
for U/t ^ 3 the system is insulating, way before one can 
invoke a 'hard' gap in the spectrum due to 'preformed 
pairs'. Fig.l.(c) highlights the low energy behaviour in 
the quasiparticle density of states (DOS). The supercon- 
ducting phase has a gap for all U and T < Tc- So does 
the large U normal state, as indicated. The low U system 
has a band like DOS for T > Tc and the most intrigu- 
ing behaviour occurs for 2 < U/t < 5, where the T > Tc 
phase has a 'dip' in the low energy DOS, i.e, a pseudogap. 

Fig.l.(d) shows how the optical spectrum varies for 
changing U and T. As the expression for cr(w) reveals, 



except at T = 0, there is in principle always low frequency 
spectral weight in the optical conductivity at any U and 
T. This arises from thermally excited quasiparticles, and 
is exponentially small when Aqp/ksT ^ 1, i.e, broadly 
in the lower right part of the U/T plane (the PI, P4 
regions). Following the same argument, PI and P3 have 
significant low frequency weight. PI and P2 of course 
also have a superfluid S function feature which P3 and 
P4 do not have. 

Resistivity: Let us shift to the resistivity. Fig. 2. From 
the expression for a(uj) it is obvious that as long as 
there is a gap in the QP spectrum, the ai term cannot 
contribute to the dc conductivity (it has a lower cutoff 
2Aqp). We tried a crude model to analyse the p{T) 
result. We assumed A^ configurations such that |Ai| is 
same at every site and is set to the T = mean field value 
appropriate to the pL and U. The phases 9i are assumed 
to be completely random and uncorrelated between sites. 
Solving the BdG equations that results from these config- 
urations, and introducing T only as in the Fermi factors, 
leads to a result that is remarkably similar to Fig. 2 for 
T > 0.3. This suggests that the detailed T dependent 
distributions of |Ai| and the spatial correlations in 9i 
are not essential for a first understanding. The p{T) for 
[/ > 6 is mainly controlled by QP activation across the 
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FIG. 3: Colour online: Temperature dependence of the QP 
DOS, N{lij) at different couplings. Panels (a)-(c) have the 
same legends, (a). Weak coupling, U/t = 2, (b). intermediate 
coupling, U/t — 6, and (c). strong coupling, U/t — 10. The 
oscillations in the DOS in panel (a) are finite size artifacts 
(even on a 32 x 32 lattice). At U/t — 2 the gap essentially 
vanishes at T ~ Tc, while ai U/t — 6 & small 'hard gap' per- 
sists to Tc and above, although lorentzian broadening gives 
the impression of a pseudogap at the highest T. For U/t — 10 
a 'hard gap' persists to T ~ 0.5 although with a clear re- 
duction with increasing temperature. Panel (d). shows the 
variation in the gap, normalised by its T = value, for the 
three couplings. The gap is inferred directly from the eigen- 
value spectrum while the DOS has a lorentzian broadening. 
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FIG. 4: Colour online: Temperature dependence of the opti- 
cal conductivity, (j{uj) in different coupling regimes, (a). For 
'weak' coupling: U/t — 2, the low oj region fills up rapidly and 
a-{u}) is almost T independent for T > Tc. (b). For interme- 
diate coupling, U/t — 6, the suppressed low cj weight persists 
to T > Tc, and the overall weight in the T > Tc spectra is 
moderately T dependent, (c). For U/t — 10 the suppression 
of low OJ weight persists to T ~ 0.5 and the integrated weight 
seems to increase slowly for T > Tc- 



T > Tc gap, while for J7 ^ 6 we observed a pseudogap 
in the spectrum generated by the toy {A^} and the p{T) 
behaviour is affected mainly by the scattering effect due 
to the random 9. The M-I transition, with growing U 
at r > Tc, is mainly due to the scattering induced by 
pairing (angular) disorder, rather than the opening of a 
clean gap in the QP spectrum. 

Density of states: The DOS has been investigated ear- 
lier too, and we confirm the known trends albeit with 
much higher resolution. Panels (a)-(c) in Fig. 3 show the 
DOS at [/ = 2, 6, 10, respectively to the left of peak, 
peak Tc, and right of peak, on the Tc{U) curve. At T = 
all plots have the usual coherence peak at the gap edges. 
These diminish rapidly with T and the gap begins to fill 
up. The fate of the gap is shown in Fig.3.(d), where 
[/ = 6 and 10 have a gap for T > (and the U = 10 
gap is strongly T dependent) while the U = 2 gap closes 



at T - Tc. This confirms that (7 ~ 3, where the M-I 
crossover occurs, will not have a clean gap for T > Tc- 

Optics: Fig. 4 shows the optical conductivity a.t U = 
2, 6, 10. We have not shown the 6 function at T < Tc 
for clarity (at T = 0, within this calculation, the entire 
weight would be in the d function). The response at 
U = 2 does not have any finite uj peak but U — 6 and 
U — 10 have a 2Aqp peak, arising from (Ti(aj), that show 
up at T ^ O.bTc, and grow with increasing T. The filling 
up of the low UJ part is however due to thermally excited 
QP's contributing via fT2(w). For T > Tc the optical 
spectral weight weakly T dependent. 

In conclusion, we have studied the attractive 2D Hub- 
bard model via a new Monte Carlo technique on large 
lattices. Our results on Tc and density of states corre- 
spond to available QMC results, but our access to real 
frequency conductivity data allows the first determina- 
tion of the metal-insulator transition in the normal state, 
and its relation to the single particle spectrum. 
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